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ABSTRACT 

The effect of a possible hydrodynamical wind on the nature of hot accretion dislc solutions is studied. 
It is found that the advection dominated branch, in the presence of a wind, maintains the self-similar 
solution for the disk structure with the temperature, 9 oc 1/r, optical depth r oc r^~^/^ and accretion 
rate m cx r^. Based on global solutions cooling due to wind energy loss and advection are found to 
be equally important. For a wide range of viscosity and wind parameters, the temperature is about 
one-tenth of the virial value and P ~ 0.9, independent of the mass accretion rate and radius. 

In the context of cooling by unsaturated Comptonization of soft photons, solutions also exist where 
radiative cooling, advection, and wind cooling are important. In this case, wind regulated solutions are 
possible. Here, the radial dependence of the critical mass accretion rate above which solutions do not 
exist is unchanged from those solutions without winds. 

The wind/advection dominated solutions are locally unstable to a new type of instability called "wind 
driven insbility" where the presence of a wind causes the disk to be unstable to long wavelength pertur- 
bations of the surface density. The growth rate of this instability is inversely proportional to the ratio of 
the radiative cooling to the gravitational energy dissipation rates and it can grow on a timescale much 
longer than the viscous timescale in the disk for sufficiently small radiative cooling efficiencies. 

Subject headings: accretion, accretion disks-black hole physics -hydrodynamics 
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1. INTRODUCTION 

Hot accretion disks are a natural site where externally 
or internally produced soft photons can be Comptonized 
to higher energies by disk electrons. It is this property 
that makes hot accretion disk models attractive in ex- 
plaining the spectra of compact X-ray binaries and ac- 
tive galactic nuclei. Alternative models to this picture 
have also been proposed involving shocks (Chakrabarti & 
Titarchuk 1995), patchy or uniform coronae (Haardt & 
Maraschi 1993; Liang & Price 1977), or spatially extended 
regions in a transition disk model (Misra, Chitinis & Melia 
1998). Detailed comparisons of the observed X-ray spec- 
tra with predictions of these different models can be used 
to determine the applicability of each model to a given 
system. In addition, constraints imposed by X-ray timing 
analysis can also be used to further distinguish between 
these models (Nowak et al 1999). Apart from these ob- 
servational tests, theoretical studies of the stability and 
self-consistencies of these models are important for deter- 
mining the circumstances under which they are physically 
realizable. 

In a seminal study, Shapiro, Lightman & Eardley (1976) 
proposed a hot accretion disk model to explain the X-ray 
spectra of the black hole system Cygnus X-1. Here, a hot 
inner disk was assumed to be cooled by external soft pho- 
tons originating from an outer cold disk. In this model, the 
gravitational energy dissipated in the hot disk heats the 
protons, and due to the inefficiency of the Coulomb inter- 
action rate in these environments, leads to a large temper- 
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ature differences between the protons and electrons. The 
outer cold disk was assumed to be transformed into a hot 
inner disk for radii less than a critical value due to a ther- 
mal instability associated with the dominance of radiation 
pressure (Lightman & Eardley 1974). Despite the suc- 
cess of this model in explaining the X-ray observations of 
Cygnus X-1, there remains several outstanding theoretical 
issues. Specifically, it was not clear whether an unstable 
cold disk would in reality evolve to a hot disk (Chen & 
Taam 1994) or whether the narrow transition between the 
hot and cold disk was energetically possible. Moreover, the 
hot accretion disk solution, itself, was thermally unstable. 

Against this backdrop, the importance of radial advec- 
tion of energy in such disks was noted by Ichimaru (1977). 
Taking this effect into account, Narayan & Yi (1994) con- 
structed a second type of hot disk solution with a self- 
similar structure called Advection Dominated Accretion 
Flows (ADAF). These solutions were, in general, radia- 
tively inefficient (with most of the energy advected into 
a black hole) compared to the non-advective, radiatively 
cooled solutions found earlier. In contrast to the radia- 
tively cooled solutions, ADAFs are thermally stable, and 
it was shown that a sharp transition between the hot in- 
ner regions and the cool outer regions could exist, provided 
an additional mechanism (e.g. conduction/turbulent dif- 
fusion) is assumed which transfers energy from the hot 
disk to the transition radius (Honma 1996). The topology 
of these solutions was studied by Chen et al. (1995) and 
Zdziarski (1998), and it was found that the advection dom- 
inated and radiatively cooled solutions only existed below 
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a critical mass accretion rate. Observational support for 
the radiatively inefficient ADAF solutions was provided by 
the X-ray spectral fitting of low luminosity galactic centers 
and quiescent X-ray Novae (Narayan et al. 1998). 

It was noted by Narayan & Yi (1994) that steady state 
rotating solutions do not exist for an adiabatic index of 
the accreting flow, Tn = 5/3. In this case the ADAF solu- 
tions become self-similar spherical Bondi accretion flows. 
This is not a general result, however, since Ogilvie (1999) 
has shown that time dependent self similar solutions for 
= 5/3 do exist. If we consider only time independent 
flows, Narayan & Yi (1994) constructed self-similar rotat- 
ing solutions in the case of a disk with a nearly equiparti- 
tion magnetic field (corresponding to an effective adiabatic 
index of 4/3). Since it is not clear that all accretion flows 
will contain such fields (e.g., Balbus & Hawley 1998), the 
applicability of these solutions could be limited. Moreover, 
Bisnovatyi-Kogan & Lovelace (2000) have recently argued 
that, in the presence of an equipartition field, the electrons 
in the disk will be preferentially heated compared to the 
protons and a two temperature disk may not exist. In con- 
trast, Blandford & Begelman (1999) have recently shown 
that a rotating self-similar flow with T]j = 5/3 is possible, 
in the presence of a wind, if the accretion flow rate in the 
disk varies with radius as M oc R^. A strong hydrody- 
namic wind from hot disks is possible, in principle, since 
the ion temperature in the disk is close to the virial value as 
shown in the early work by Takahara, Rosner, & Kusunose 
(1989) and Kusunose (1991). Recently confirmation of 
these results has been obtained by Beckert (2000) for a 
variety of viscosity prescriptions. In a parallel develop- 
ment, Chakrabarti (1999) and Das (1999) have studied the 
possibility of outflows in the context of shock/centrifugal 
barrier models for the accreting flow. Comparisons with 
observations reveal that the X-ray spectra of such a wind 
driven self-similar flow called ADIOS (Advcction domi- 
nated inflow-outflow solutions) can explain the observed 
spectra of black hole candidates in quiescence (Quatacrt 
& Narayan 1999). In addition, Lin, Misra, & Taam (2000) 
find that the spectral characteristics of high luminosity 
black hole systems suggest that winds may be important 
for these systems as well. 

Theoretically, the conditions under which wind dom- 
inated accreting fiows exist are unclear. For example, 
Abramowicz, Lasota, & Igumenshchev (2000) argue, based 
on the sign of the Bernoulli function and on detailed two 
dimensional hydrodynamical simulations that winds are 
not expected for disks characterized by low viscosities 
{a < 0.1). Although modest outflows have been seen 
in simulations of disks characterized by high viscosities 
(Igumenshchev & Abramowicz 1999; see also Stone et al. 
1999), the relaxation of the assumption of identical adi- 
abatic indices in the disk and wind (reflecting different 
mechanisms for driving accretion and outflow) can en- 
hance the effectiveness of winds (see Chakrabarti 1999). 
The existence of winds can also be affected by the pres- 
ence of convection which can modify the temperature pro- 
file (Narayan, Igumenshchev, & Abramowicz 2000). In 
particular, they find that the thermal disk structure is de- 
pendent on the magnitude of the a viscosity parameter and 
the direction in which angular momentum is transported 
by the convection. For outward transport of angular mo- 



mentum the ADAF solutions are recovered, whereas for 
inward transport both non-accreting solutions with low- 
ered temperatures {a < 0.05) and the standard ADAF so- 
lutions are found {a > 0.05). Thus the effect of convection 
on the existence of winds remains inconclusive especially 
since it is dependent on the unknown distribution of en- 
ergy dissipation in the vertical extent. It is clear that a 
strong wind could exist under certain conditions, and it 
would be useful to study the effects of such a wind on the 
structure of hot accretion disks. 

In ADIOS solutions it is generally assumed that the 
mass flow rate has a power-law dependence on radius with 
the power law index, P, treated as a parameter (Blandford 
& Begelman 1999; Beckert 2000). In principle, the vari- 
ation of the accretion rate with radii should be obtained 
self-consistently from the surface wind velocity and the 
specific energy of the wind matter. This, in turn, depends 
on the disk structure, geometry and the thermal proper- 
ties of the wind. A realistic computation of this disk/wind 
configuration is difficult due to the complexity of the prob- 
lem and the unknown vertical distribution of the viscous 
dissipation in the disk. However, we show that significant 
progress can be attained if only the dependence of the 
surface velocity on the disk structure is estimated and its 
numerical value is parameterized to take into account the 
uncertainties involved. 

In the next section, the basic equations of hot disks with 
winds are described. In §3 the results of the calculations 
are presented and in §4 the stability of the solutions is dis- 
cussed. Finally, we summarize the results and conclude in 
the last section. 

2. DISK STRUCTURE EQUATIONS 

The standard equations for a Newtonian, keplerian hot 
accretion disk are modified to take account of mass, en- 
ergy, and angular momentum loss due to a hydrodynami- 
cal wind. To facilitate the presentation of these equations 
in a compact form, the following dimensionless quantities 
are introduced: the ion temperature 9 = kTi/mpC^; the 
electron temperature Og = kTe/mgC?; the radius r = R/rg 
where = GMjc? and M is the mass of the black hole; 
the accretion rate m = M/MEdd where the Eddington ac- 
cretion rate is defined to be Msdd = ^'KrUpCrg/aT; Thom- 
son optical depth r = naxH, where n and H are the 
number density and half-height of the disk respectively; 
the speed of the wind at the disk surface Vw = Vw/cg 
where Cg = cO^/"^ is the local isothermal sound speed 
of the disk and the Coulomb energy exchange frequency 
V* = {rglc)Vc. 

For a non-magnetized two-temperature disk i.e. when 
the ion temperature (Tj) is much higher than the electron 
temperature (Te), the local pressure is simply Pg = nkTi 
and the equation for hydrostatic equilibrium can be writ- 
ten as 

{H/Rf = 9r (1) 
In the presence of a wind the accretion rate is no longer 
constant and hence mass conservation gives, 

^ = ATTRmpuV^ (2) 
art 

which in dimensionless form can be re- written using equ. 
(1) as. 
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The torque equation (e.g. Bcckcrt 2000) is, 

where T is the viscous torque and where it has been as- 
sumed that the specific angular momentum of the wind is 
same as that in the disk i.e. I = \/GMR. The viscous inte- 
grated stress is described in terms of the a-prescription of 
Shakura & Sunyaev (1973), and it scales as the pressure, 
i.e. T = aPgH for a disk in hydrostatic equilibrium. The 
above equation can then be written in dimensionless form 
as 

d(j^ = lrn(r)r-V^ (5) 

Note that since the accretion rate has a radial dependence 
(duo to the presence of a wind) the above differential equa- 
tion cannot be trivially integrated. 

For a disk in thermal equilibrium, the local vertically in- 
tegrated energy dissipation rate {Qd) should be equal to 
the sum of the Coulomb energy loss {Qc), the advective 
energy flux (Qa) which is non-local radial energy transfer 
rate associated with the flow and the energy flux carried 
away by the wind {Qw)- The energy dissipation rate due 
to the viscous torque for each half of the disk is 



(6) 



A fraction of this dissipated energy will be transfered by 
the protons to electrons via Coulomb interactions. The 
Coulomb energy exchange flux is 

3 MEddC^ me 

QC = -^ 15 '^^ci^ ^e) (7) 



3 MEddC^ 



rrij, 



(8) 



Here, the dimensionless Coulomb frequency is given by, 

U* = ^U, = 4.4 X 10-4lnAT^-3/2^-V2r-3/2 (g) 

where InA « 15 is the Coulomb logarithm (Spitzer 1962). 
The advective energy flux is 

M Pg 1 MEddC^ . nc -2 
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where the advective factor is defined to be (Chen et. 
al. 1995) 

3r dO r dr r dH 

It has been implicitly assumed in this definition that the 
specific energy of the gas is {3/2)kTi or in other words, the 
adiabatic index of the accretion flow Td = 5/3. Note that 
self-similar (ADAF) solutions without winds, oc r~^, 
H (X r and r oc r~^/^, leads to ^a = 0. This is consis- 
tent and expected since, for Td = 5/3 and in the absence 
of a wind, there are no steady self-similar (ADAF) rotat- 
ing solutions and only spherical Bondi accretion flows are 
allowed (Narayan and Yi 1994). 

The energy flux carried away by the wind is 



Qy^^ = nVwEw = z 5 — TVu,ewr °' 

Sir ri 



(12) 



Here Eyy is the average energy imparted to a proton in the 
wind from the local disk and Cu, = Ew/{GMmp/2R) is 
its dimensionless form. The external energy required for a 
particle to escape to infinity is w GMmp/2R~ (3/2)(fcTi — 
kTf) where Ti and Tf arc the proton temperature of the 
wind at the disk surface and at infinity respectively. If the 
wind is isothermal or if fcTj << GMmp/2R then this ex- 
ternal energy input per particle is simply GM/2R or the 
binding energy of the particle. Further, if this energy is 
obtained only from the local region of the disk from where 
the particle originates, then ~ 1. Thus, < 1 corre- 
sponds to the case where only a fraction of the required 
energy is obtained from the local region of the disk and 
the rest presumably is due to heating of the wind from a 
different region of the disk. The > corresponds to 
the unlikely case where more energy than required is de- 
posited from the local disk to the wind. This additional 
energy may be radiated away by a radiativcly cooled wind 
and/or appear as kinetic energy of the wind particles at 
infinity. The uncertainties involved in the energy coupling 
between the wind and the disk is factored into which 
is treated as a parameter in this work. 

The equation describing conservation of energy [Qd = 
Qc + Qa + Qw ) can be written in a compact form using 
equs. (6), (8), (10) and (12) as 

u*r^/'^ ^ 2 mjA ^ 1 Vwe^ _ ^ 
a 3 arr^/^ 3 a9r 
The value of the electron temperature (^e) which is re- 
quired to calculate the Coulomb frequency (equ. 9), de- 
pends on the dominant radiative cooling mechanism op- 
erating in the hot disk. Following Zdziarski (1998) it is 
assumed here that the disk cools primarily by unsaturated 
Comptonization of soft photons. This is motivated by the 
high energy spectral characteristics of black hole binaries 
in the low-state since the Comptonization process is char- 
acterized by the Compton y-parameter {y = AkTe/meC^) 
max(r, r^) w 1. Thus 

4maa;(r, r^) ^^^^ 
can be used to estimate the electron temperature using 
y as a parameter of order unity. Technically the electron 
temperature profile in the disk should be obtained by solv- 
ing the geometry dependent radiative transfer equations in 
both vertical and radial directions. Thus it should be em- 
phasized that equation (14) is a simplifying assumption 
which may be approximately valid only under certain con- 
ditions (Shapiro, Lightman & Eardley 1976). 

Determining the wind velocity at the surface of the disk 
is complicated because it depends on the vertical structure 
of the disk which in turn depends on the detailed varia- 
tion of the unknown viscosity with height. Furthermore, 
the opening angle for the wind flow needs to be determined 
using the geometry of the accreting flow and the pressure 
gradient in the wind in the radial direction. However, a 
qualitative estimate of the dependence of the wind velocity 
with respect to the average temperature and gravitational 
potential in the disk, can be obtained by studying the 
simpler geometry of a spherical flow from a radius R, tem- 
perature Ti and potential « GM/2R. The factor of two 
in the potential takes into account the angular momentum 
of the disk. The estimate obtained from such a spherical 
flow can be scaled by an unknown factor r] which would 
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then take into account the uncertainties in the geometry 
and vertical structure of the disk, 

Vyj='qvs{0r) (15) 

where Vg is the velocity of the wind at the surface of the 
sphere divided by the isothermal sound speed (= 
Thus the motivation here is to obtain a qualitative depen- 
dence of the disk surface wind velocity with temperature 
and radius of the disk, rather than a quantitative esti- 
mation of its magnitude. The standard equations of the 
hydrodynamical theory of stellar winds have been used to 
compute the wind velocity at the surface of the sphere. Us- 
ing Bernoulli's equation and the fact that the Mach num- 
ber at the sonic point is unity, the following relation for 
the sonic radius is obtained 



r 

r-i 



A 



(16) 



where Tc is the location of the sonic point, is the ra- 
dius of the sphere, T is the adiabatic index of the wind 
and A is the ratio of the gravitational potential at the sur- 
face to the temperature of the sphere i.e. A = <p/kTi = 
GMnrip/ {2rskTs). Note that the adiabatic index in the 
wind may differ from that in the disk (see Chakrabarti 
1999). Conservation of mass and the equation of state, 
Pg oc (where p is the mass density) gives. 



.A 1 r+i 
= (-)2 r-i 



(17) 



In this work, equations (16) and (17) are solved simulta- 
neously to yield Ws(A, F) where A for an isothermal disk is 
1/(26*0. If « A, then 



Vs{\T) 
where Fi = 



U2[(5-3F)(1 



A(r-i) 
r 



)]^^(18) 



3r-5 1 p _ 7r-9 
ly and i 2 - 2(FrT) 



2(r 



For the isothermal 



case (F — > 1), equ. (18) reduces to 

t;,(A,F^l) = (^)2exp(3/2-A) 



(19) 



There are several points which should be noted here. First, 
equs. (18) and (19) seem to suggest that Vs has a max- 
ima for variations with A. This is an artifact of the as- 
sumption t;^ << A and if equs. (16) and (17) are used 
instead (as is done in this work), then Vg is a monotonic 
function of A. Second, there is a maximum value of A 
beyond which there arc no wind solutions or = for 

< Omax = (r - l)/(2Fr). Third, inspection of equ. (18) 
reveals that for 9 <• 9„iax , is a rapidly varying function 

01 6, a factor two increase in 9 could lead to an increase of 
two orders of magnitude in Vw (Lin, Misra, & Taam 2000). 

The equations corresponding to conservation of mass 
(3), angular momentum (5) and energy (13) can be solved 
to obtain the disk structure using the estimated values for 
electron temperature (14) and wind velocity (15). 

3. DISK SOLUTIONS WITH WINDS 

3.1. Local Solutions 

The differential equations governing the disk structure 
can be solved using appropriate boundary conditions. 
However, these equations can be reduced to algebraic ones, 
by treating the logarithmic derivatives of physical quan- 
tities as constant (i.e., independent of radii) parameters. 



The solutions for this set of algebraic equations depends 
only on the local radius and hence are local solutions. Such 
solutions generally provide important insight on the disk 
structure and its dependence on disk parameters. 

To obtain local solutions it is assumed that the logarith- 
mic derivative of temperature (^t = —dlogT/dlogR) and 
accretion rate (P = dlogM/dlogi?) with respect to the 
logarithm of the radius are constants. Integrating equa- 
tion (5) with the constant of integration chosen such the 
viscous torque vanishes in the inner radius, leads to 



m(r)r 2j(r-) 
(1 + 2P)a9 



(20) 



where J (r) = 1 — (6/r)^+^/^ takes into account the torque 
free inner boundary condition at r = 6. The advective 
fraction (equ. 11) is reduced to 

u = mT-i)+p (21) 

where it has been assumed that J(r) « 1 is a constant. 
While it is convenient to consider as a parameter, P can 
be obtained self-consistently by mass conservation (equ. 
3), 

.1/2 



P = 



TV 



m 



-r]Vs{9r) = 



1 



ri.,Vs{9r) 



l + 2P'-a 
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9r 



)• 



(22) 



Hence, the local disk solution is obtained by solving the 
energy equ. (13) with equs. (20), (21) and (22). Figure 
1 shows the variation of the normalized accretion rate, m, 
with optical depth, t, for r = 40. Comparison with Fig- 
ure 1 of Zdziarski (1998) shows that the topology of local 
solutions does not change with the inclusion of a wind. 
For accretion rates less than a critical value (m < merit), 
there exists a radiatively cooled solution branch and an 
advection/wind dominated solution branch. No solutions 
exist for rh > rticrit- Note, that for F = 5/3, comparisons 
cannot be made with respect to solutions without winds 
since steady state solutions do not exist. 

For the advection/wind dominated solution branch, the 
Coulomb exchange rate can be neglected {Qc ~ 0) and 
equ. (13) can be written for = 1 as 



|=P(l + 2P)^r+f (^)(^) 
2 2 a ur 

= {^)vs{9r){l + ^) 
a 29r 



(23) 



The above with equs. (20) and (22) lead to self-similar 
solutions for an advection/wind dominated flow: 9 oc 1/r, 
r oc r^~^/'^ and rh oc r^. Thus, with the inclusion of a 
wind, self-similar flows are maintained which is a direct 
consequence of the form of the wind velocity chosen in 
equ. (15). For example, if the parameter 77 is a function 
of radius rather than a constant as assumed here, self- 
similar solutions will no longer be allowed. The solution 
to equ (23) is shown in Figure 2, where 9r is plotted as a 
function of r^/a for different values of F and Ew . For val- 
ues oi T]/a > 2, 9r depends only logarithmically with r]/a 
and is ~ 0.1. The ratio of the wind to advective cooling 
Qw/Qa = sw / {'29r) w 5 making both cooling channels 
equally important for this solution branch. It is illustra- 
tive to solve for P using equs. (22) and (23) yielding 



P = ( - 1 + Vl + 24/(eM. + 2^r))/4 



(24) 
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which impUcs that P « 0.9 for ew ~ 1 and 9r Ri OA. This 
is shown in Figure 3 where P is pfotted as a function of 
rj/a. Note that from equ. (24) a smaller value of ew = 0.5 
corresponds to P « 1.3 and a larger value ew = 1-5 corre- 
sponds to P w 0.7. Thus P remains of order unity as long 
as ew ~ 1 which is expected unless the wind is predom- 
inately heated non-locally or if the energy transferred by 
the disk to the wind greatly exceeds the energy required 
to generate the wind. 

At the critical accretion rate, rncru, radiative cooling 
becomes comparable to the local dissipation (Qc Qd)- 
Using the self-similar solution obtained above, rhcrit oc 
^i^ic}i is the same dependence obtained for hot solu- 
tions without winds (Zdziarski 1998). The value of rhcrit is 
however, generally lower than for solutions without winds. 

In summary, for a wide range of values for the vis- 
cosity parameter {a) and the wind parameter (77), ad- 
vection/wind dominated solutions arc characterized by 
9 ~ 0.1/r and P « 0.9, independent of the accretion rate 
or radius. We note that the result that the ion tempera- 
tures of wind dominated solutions are about an order of 
magnitude smaller than the no-wind ADAF solutions jus- 
tifies the assumption of a keplerian disk. 

3.2. Global Solutions 

It is prudent to confirm the inferences obtained from 

local analysis by constructing global solutions of the disk. 
These global solutions involve solving the three differential 
equations for conservation of mass, angular momentum 
and energy (equs. 3, 5 and 13) with three independent 
boundary conditions. For convenience, the boundary con- 
ditions adopted in this work are that the viscous torque 
vanishes at the inner radius (rj = 6), and the accretion 
rate, m(ro), and the temperature, 9{ro), are specified at 
the outer radius, rg- For a given rh{ro), 6{ro) is computed 
from the local solution described in §3.1. 

Figure 4 illustrates the variation of the mass accretion 
rate with respect to radius for different values of m{ro)- 
Here, 9{ro) is chosen from the advection/wind dominated 
branch of the local solutions and hence, the global solution 
presented is also advection/wind dominated. The global 
solutions are, indeed, similar to the type inferred from lo- 
cal analysis with P « 0.9 except when r ^ 6 where as 
expected the flow is no longer self-similar. The dotted line 
is the critical accretion rate (beyond which hot disk solu- 
tions no longer exist) computed using local analysis but 
with values of P and £_t taken from the global solution. 

Figure 5 presents the global solutions for the case when 
9{ro) is chosen from the radiatively cooled branch of the lo- 
cal solutions. For m(ro) << mcrit(ro), the accretion rate is 
a constant since winds are not important (i.e., 9 << 1/r). 
However, for m(ro) ^ fhcritiro), winds do become im- 
portant for certain regions of the disk. In such cases the 
wind regulates the flow such that m(r) ^ rticritir)- These 
wind regulated solutions occur for a limited range of m{ro) 
which is due to the weak dependence of rhcrit with radius 
[merit oc r^/^^). In contrast, the accretion rate at the in- 
nermost regions of the disk can be significantly reduced 
from m(ro) for the advection dominated branch. For ra- 
diative cooling mechanisms other than soft photon Comp- 
tonization, rhcrit may depend more steeply on radius, lead- 
ing to a wider range of situations where a wind-regulated 



accretion flow would exist. 

4. STABILITY OF DISK SOLUTIONS WITH WINDS 

Radiatively cooled hot disks (without winds) are known 

to be thermally unstable. This result remains valid for 
the radiatively cooled branch described in §3.1, simply be- 
cause winds are not important for these solutions. Advec- 
tion dominated hot disks (without winds) are known to be 
thermally stable and again this result is not altered for the 
advection/wind dominated solution since the energy loss 
rate due to the wind increases with temperature. 

Formally the criterion for thermal stability is that the 
derivative of the energy dissipation rate with respect to 
temperature should be smaller than the derivative cooling 
terms with respect to temperature, at a constant surface 
density. This leads to (Appendix: equ. Al) 



St = + 2(5a + S,ygC 



ID 



> 



(25) 



as the criterion for thermal stability provided that the 
advection timescale is longer than the thermal timescale. 

Here, = dlogVs/d\og9 > 1. For the wind/advection 
dominated solutions, Qd = Qa + Qw and Qc << Qd, 
hence the disk will be locally stable to thermal perturba- 
tions. For the radiatively cooled branch, at high accretion 
rates, advection and wind cooling may be sufficiently im- 
portant to stabilize the disk. 

On timescales longer than the viscous timescale, the sur- 
face density (or optical depth) will respond to fluctuations 
in the disk and its temporal behavior is governed by, 

where Vr is the radial velocity. From the torque equation 
(see equ. 5), it can be seen that fluctuations in Vr are 
linearly proportional to variations in temperature in the 
long wavelength limit. On the other hand Vw/H{(x Vw) 
is a rapidly increasing function of temperature. Thus, in 
thermal equilibrium, if r is inversely proportional to tem- 
perature, the disk will be unstable to long wavelength per- 
turbations. In other words, if a decrease in r leads to an 
increase in the temperature, the mass loss rate by the wind 
will increase rapidly which will cause a further decrease in 
T leading to a runaway situation. 

Technically, for the local wind dominated solution de- 
scribed in §3.1 without radiative cooling (Qc = 0), the 
temperature of the disk is independent of r (equ 23). How- 
ever, even if there is only a small contribution due to radia- 
tive cooling (i.e. Qc << Qd), an increase in t will cause 
Qc to increase which leads to a decrease in the tempera- 
ture. Thus, wind dominated solutions are subject to the 
instability described above in the long wavelength limit. 
A linear stability analysis leads to the dispersion relation 
for A >> P (Appendix: equ. A4) 



5 T4.N/QcND/t i\ 



(27) 



2'P''S't 

For wind dominated solutions, S^yg > 1 and the disk is 
thermally stable {St > 0), hence iiv is real and positive 
indicating that the mode is unstable. For a radiatively 
cooled disk, St < and this mode is stable for such so- 
lutions. The growth time of the wind instability (« ^/i^) 
can be significantly longer than the viscous (or accretion) 
timescale (« R/Vr) if Qc/St « Qc/Qd « 1- This im- 
plies that the disk structure can be described in terms of 
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local steady state solutions during the linear growth of the 
wind driven instability. Equ. (27) and the inferred argu- 
ments depend on the form of the radiative cooling mech- 
anism (equ. 14) assumed in the analysis. In the general 
case, Qc oc , and it can be shown that the disk will suf- 
fer this instability provided /? > 1 and if St > 0. For the 
unsaturated Comptonization model (i.e. equ. 14) /3 = 7/2 
and hence the disk is found to be unstable. Note however, 
that /3 < 1 corresponds to the situation where the elec- 
tron temperature is proportional to the optical depth and 
varies faster than 6^ oc t^/^ . Since for most physical radia- 
tive processes, Og is expected to be inversely proportional 
to T, hot disks with winds would in general be unstable. 
The above results are based on local and linear stability 
analysis and require confirmation by global simulations. 

5. SUMMARY AND DISCUSSIONS 

The effect of a hydrodynamical wind on a hot accre- 
tion disk has been studied for disks characterized by an 
adiabatic index, Td — 5/3. It is foimd that the topol- 
ogy of local ADAF solutions found earlier (for Tn < 5/3) 
remains the same for To = 5/3 with the inclusion of a 
wind. In particular, there exists two branches of solutions 
for accretion rates below a critical value, while there are 
no solutions for higher accretion rates. One solution corre- 
sponds to the radiative cooling branch where the effects of 
the wind arc negligible whereas both winds and advection 
are important for the second branch. 

We have shown that the local solutions for the 
wind/advection dominated branch may be described in a 
self-similar form with the temperature, ^ oc 1/r, optical 
depth T cx r^~^/'^ and accretion rate rh cx . For a wide 
range of viscosity and wind parameters, the temperature 
is about 10% of the virial value {6r « 0.1) and P « 0.9, in- 
dependent of the accretion rate and radius. These results 
were confirmed by global analysis for which self-similar so- 
lutions were obtained for radii much greater than the inner 
disk radius, thereby providing a theoretical basis for the 
model advanced by Blandford & Begelman (1999). Global 
solutions also revealed that under certain boundary con- 
ditions (see §3.2), the accretion rate m « merit r^/^^. 
For these "critical" solutions radiative cooling and energy 
loss due to the wind, were both important and hence they 
are not self-similar. 

The wind/advection dominated solutions arc thermally 
stable, but are locally unstable to a new type of insta- 
bility which we denote as "wind driven instability." The 
presence of a wind causes the disk to be unstable to long 
wavelength perturbations in the surface density. Since 
the growth rate of this instability is inversely proportional 
to the ratio of the radiative cooling to the gravitational 
energy dissipation rates, the timescale for the instability 
could be much larger than the viscous time scale in the 
disk for a low radiative cooling efficiency. 

The accretion rate index, P, w 1 implies that, if the 
outer radius of the hot disk is much larger than the inner 
one, most of the accreting matter is lost in the form of 
a wind rather than accreted onto the star. If the com- 
pact object is a neutron star or white dwarf, it is this 
residual accretion which gives rise to the X-ray luminosity 
from the surface of the star. Loeb, Narayan & Raymond 
(2000) point out that in such a scenario (with nearly simi- 



lar outer boundary conditions), the surface emission from 
the star will scale roughly as oc Rf~^, where i?* is the 
radius of the star. From the observed luminosities of soft 
X-ray transients and cataclysmic variables, they estimate 
P = 0.9 ± 0.5, which is in agreement with the theoretical 
results presented here. 

For this radial dependence of the mass flow rate in the 
disk, the total energy dissipated in the disk scales as log 
(Ro/Ri), where Ro and Ri are the outer and inner radii 
of the hot disk, rather than R~^- Thus both the inner 
and outer regions of hot disk are equally important ener- 
getically. It is generally assumed that the hot disk is sur- 
rounded and connected to an outer cold disk by a narrow 
transition region. The energy required to make this tran- 
sition cannot be generated locally if the transition zone 
is small (Misra & Melia 1996) and is probably due to 
global transfer from the inner hot regions to the transi- 
tion zone (Honma 1992, see also Abramowicz ct al. 2000). 
For the wind/advection dominated case this may not be 
possible since the energy generated in the inner regions 
is only comparable to that of the outer region. In this 
case, the transition zone would then have to be extended 
and could significantly contribute to the luminosity of the 
source (Misra & Melia 1996). However, the above spec- 
ulation has to be verified by detailed multi-dimensional 
simulations of the transition zone taking into account the 
possibility that the hot disk may be wind dominated. 

Since the radiative cooling rate (equ. 8) for the self- 
similar solution scales as oc ij(i4^-i5)/4 ij-i/^ for 
P « 0.9, most of the luminosity arises from the outer 
regions of the hot disk. This result, however, depends on 
the radiative cooling mechanism operating in the hot disk. 
For example, if the electron temperature is chosen to be 
9e OC instead of equ. 14, then the radiative flux is 
oc i?2P-3-3ae/2 „ j^-(i+3aj2)_ ^^^^^ > 2/3 then 

the luminosity will be dominated by the inner region. 

The terminal velocity of the wind should be approxi- 
mately equal to the sound speed at the surface of the disk. 
This implies that the bulk energy per proton of the wind 
flow at infinity is « O.IGM/R since the temperature of 
the disk 9 ^ O.lr For P « 1 most of the wind arises from 
the outermost radius and hence the total bulk energy flux 
in the wind w O.lMwindGM/ Ro, where Mwind is the total 
mass loss rate and Ro is the outer radius of the hot disk. 
Since the radiative efficiency of such disks is << 0.1, the 
bulk energy of the wind flow in some circumstances may 
be much larger than the radiative flux from the source. 

The viscous timescale for the self-similar solution is 
_R ^ r 

{l + 2P){er)a^ c> ^ 
which numerically is ty ~ 250(rg/c) ~ 0.01 s for a = 0.1, 
R = 50r(, and a 10 Mq black hole. Thus, a wind 
driven instability may occur on timescales ~ 1 — 10 s 
for Qc/Qd ~ 0.1 — 0.01. The instability may not lead 
to disruption of the disk, but instead give rise to large 
amplitude oscillations. For example, the surface density 
could increase, accompanied by a decreasing temperature. 
When the temperature decreases sufficiently, the wind is 
quenched and radiative cooling becomes important. Since 
the radiatively cooled disk will be thermally unstable, the 
system may return to the wind dominated state on thermal 
timescales, giving rise to cyclic behavior. Alternatively, 
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the surface density may decrease and the disk may then be 
evacuated and finally filled again on the viscous timescale. 
In either case, the disk may evolve through epochs of high 
wind activities on timescale w seconds. Such a disk/ wind 
interaction may be important for explaining the spectral 
and temporal behavior of the regular bursting activity ex- 
hibited by the microquasar GRS 1915-105 (Taam et al 
1997). 

Future work involving time dependent numerical simu- 
lations of hot disk with winds arc needed to confirm the 
above speculations and to quantify the rate of mass loss 
and varying spectral characteristics during such events. 



Such studies may also shed light on the relevance of wind 
instabilities as a mechanism for either oscillatory behavior 
or evaporative behavior in the hot inner disk regions. Be- 
sides these issues, further theoretical investigations of the 
transition from a cold to a hot wind dominated disk and 
the role of two dimensional effects along the lines suggested 
by Turolla & Dullemond (2000) would also be illuminat- 
ing. 

R.M. gratefully acknowledges support from the Lind- 
heimer Fellowship at Northwestern University. 
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APPENDIX 

LINEAR STABILITY ANALYSIS IN THE LONG WAVE-LENGTH APPROXIMATION 

The linear response of the disk structure is calculated for perturbation of the ion temperature 9 9{1 + S9e^^"*~'^^^) 
and optical depth r — > r(l + Jre'^'^*"'^''^), where 66 and 5t are small dimensionless quantities. In the long wavelength 
approximation, it is assumed that fcr << 1 such that e''^'" « 1 which implies that the radial derivatives of perturbed 

quantities can be neglected. 

In this formalism, the linear response of physical quantities (e.g., the energy flux) can be obtained directly from the 
dependence of the quantity on temperature and optical depth. The dissipative energy rate Q d oc 9t (equ. 6) and hence 
SQd = S6 + 5t. The Coulomb energy exchange flux Qc oc t'^I'^9^1'^ (equs. 8, 9, 14) while the energy flux carried away 
by the wind Qw cx Tv^ui9) (cqu. 12 ). The advcctivc energy flux (cqu. 10), Qa oc m9 oc t9'^ since m cc t9 (cqu. 5) for 
a self-similar steady state solution with kr « 1. Note that in the long wavelength approximation the advective factor 
{^a) is invariant to perturbations. Thus the normalized response of the total energy flux Qt = Qd — Qc — Qw — Qa is, 

1 7 

QtSQt = {Qd - ^Qc - ^voQw - 2Qa)S9 + {Qd - ^Qc - Qw - Qa)St 

= -StS9 - ^QcSt (A1) 

where = d\ogVs/dlog6. On thermal timescales 5t « 59 and —QtSQt/S9 = S't > is the criterion for thermal 
stability (equ. 25). 

The time dependent form of the mass conservation equation (3) governs the temporal behavior of the surface density 
or optical depth (equ. 26) 

.Tasdr 1 dm 0/9 , , „, 

-^7 = -^-^^-^ ^2 
c at r or 

Noting that ifi oc t9, the linearized form of the above equation is 

i^r5T={-){-'^)59{l-U) (A3) 
rg r or 

For timescales much longer than the thermal timescale, 6Qt ~ i.e. the system will be in thermal equilibrium. Combining 
equations (Al) and (A3), yields the dispersion relation in the long wavelength limit (equ. 27) 

^'-=^(|)(^)^'(C..-1) (A4) 
where Vr is the magnitude of the radial velocity at radius R. 
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Fig. 1. — Variation of the normalized accretion rate (rri) with optical depth (t) for r = 40. Solid line: a = 0.5; Dashed line: a = 1.0. Other 
parameters are r] = I, ey/ = 1, j; = 1, F = 1.2 and = 1. 
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Fig. 2. — Variation of dr with rj/a for F = 1.01, ew = 1 (solid line), F = 1.01, ew = 0.5 (short dashed line) and F = 1.2, ew = 1 (long 
dashed line). 
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Fig. 3. — Variation of the accretion slope P with rj/a for F = 1.01, ew = 1 (solid line), F = 1.01, ew = 0.5 (short dashed line) and 
F = 1.2, ew = 1 (long dashed line). 
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Fig. 4. — Accretion rate vs. radii for different values of accretion rate at the outer boundary (r = 400). The temperature at the outer 
boundary is chosen from the advection/wind dominated branch of the local solutions. Dashed curve represents the local critical accretion 
rate. Parameters used for the plot are a = 0.1, y = 1, rj = I, ew = 1 and T = 1.2 




Fig. 5. — Accretion rate vs. radii for different values of accretion rate at the outer boundary (r = 400). The temperature at the outer 
boundary is chosen from the radiatively cooled branch of the local solutions. Dashed curve represents the local critical accretion rate. 
Parameters used for the plot are a = 0.1, y = 1, rj = 1, ew = 1 and T = 1.2 



